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The lattice approach to QCD facilitates the numerical evaluation of expectation 
values without recourse to perturbative techniques. Although the lattice formula- 
tion is almost as old as QCD itself and first simulations of the path integral have 
been performed as early as in 1979 [0], only recently computers have become pow- 
erful enough to allow for a determination of the infinite volume continuum light 
hadron spectrum in the quenched approximation to QCD within uncertainties of 
a few per cent. To this accuracy the quenched spectrum deviates from experi- 
ment. Some collaborations have started to systematically explore QCD with two 
flavours of light sea quarks and the first precision results indeed indicate such 
differences. 

Lattice QCD is a first principles approach; no parameters apart from those 
that are inherent to QCD, i.e. a strong coupling constant at a certain scale and 
n quark masses, have to be introduced. In order to fit these n + 1 parameters, 
n + 1 low energy quantities are matched to their experimental values: the lattice 
spacing a(g,rrii), that results from given values of the bare lattice coupling g 
and (in un-quenched QCD) quark masses m^, can be obtained by fixing m p as 
determined on the Lattice to the experimental value. The lattice parameters 
that correspond to physical m u ~ rrid can then be obtained by adjusting m n /m p ; 
the right m s can be reproduced by adjusting mx/^ip or m^/m p to experiment 

*Plenary talk presented at "Nuclear and Particle physics with CEBAF at Jefferson Lab", 
Dubrovnik, November 3-10, 1998. 
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etc.. Once the scale and quark masses have been set, everything else becomes 
a prediction. Due to the evaluation of path integrals by use of a stochastic 
process, lattice predictions carry statistical errors which can in principle be made 
arbitrarily small by increasing the statistics, i.e. the amount of computer time 
spent. In this sense, it is an exact approach. 

Lattice results have in general to be extrapolated to the (continuum) limit 
a —>■ at fixed physical volume. The functional form of this extrapolation is 
theoretically well understood and under control. This claim is substantiated by 
the fact that simulations with different lattice discretisations of the continuum 
QCD action yield compatible results after the continuum extrapolation. For 
high energies, an overlap between certain quenched Lattice computations and 
perturbative QCD has been confirmed too [0], excluding the possibility of fixed 
points of the /5-function at finite values of the coupling, other than g = 0. After 
taking the continuum limit an infinite volume extrapolation is performed. Results 
on hadron masses from quenched evaluations on lattices with spatial extent La > 
2 fm are virtually indistinguishable from the infinite volume limit within typical 
statistical errors in most cases. However, for QCD with sea quarks the available 
information is not yet sufficient for definite conclusions, in particular as one might 
expect a substantial dependence of the on-set of finite size effects on the sea quark 
mass(es). 

The effective infinite volume limit of realistically light pions cannot be re- 
alised at a reasonable computational cost, neither in quenched nor in full QCD. 
Therefore, in practice another extrapolation, in the quark mass, is required. This 
extrapolation to the correct light quark mass limit is theoretically less well under 
control than those to the continuum and infinite volume limits. The parametri- 
sations used are in general motivated by chiral perturbation theory and the the- 
oretical uncertainties are the dominant source of error in latest state-of-the-art 
spectrum calculations. 

Many important questions are posed in low-energy QCD: is the same set 
of fundamental parameters (QCD coupling and quark masses) that describes 
for instance the hadron spectrum consistent with high energy QCD or is there 
place for new physics? Are all hadronic states correctly classified by the naive 
quark model or do glueballs, hybrid states and molecules play a role? At what 
temperatures/densities does the transition to a quark-gluon plasma occur? What 
are the experimental signatures of quark-gluon matter? Can we solve nuclear 
physics on the quark and gluon level? Clearly, complex systems like iron nuclei 
are unlikely ever to be solved from first principles alone but modelling and certain 
approximations will always be required. 

It is desirable to test model assumptions, to gain control over approximations 
and, eventually, to derive low-energy effective Lagrangians from QCD. Lattice 
simulations are a very promising tool for this purpose and in the first part of this 
article I will try to give a flavour of such more theoretically motivated studies 
before reviewing recent results on glueballs and exotic hybrid mesons as well as 
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discussing the light hadron spectrum. 



2 The confinement scenario 

Two prominent features of QCD, confinement of colour sources and spontaneous 
breaking of chiral symmetry, are both lacking a proof. They appear to be related, 
however: in the low temperature phase of broken chiral symmetry, colour sources 
are effectively confined. Clearly, an understanding of what is going on should help 
us in developing the methods required to tackle a huge class of non-perturbative 
problems. 

It is worthwhile to consider the simpler pure SU(N) gauge theory. In this 
case confinement can be rigorously defined since the Polyakov line is an order 
parameter for the de-confinement phase transition that is related to spontaneously 
breaking a global Zn symmetry. I will present some results that have been 
obtained in the computationally cheaper SU(2) gauge theory whose spectrum 
shares most qualitative features with that of SU(3). 

In the past decades, many explanations of the confinement mechanism have 
been proposed, most of which share the feature that topological excitations of 
the vacuum play a major role. These pictures include, among others, the dual 
superconductor picture of confinement || |J and the centre vortex model ||. 
Depending on the underlying scenario, the excitations giving rise to confinement 
are thought to be magnetic monopoles, instantons, dyons, centre vortices, etc.. 
Different ideas are not necessarily exclusive. For instance, all mentioned excita- 
tions are found to be correlated with each other in numerical as well as in some 
analytical studies, such that at present it seems to be rather a matter of personal 
preference which one to consider as more fundamental. 

Recently, the centre vortex model has enjoyed renewed attention ||. In this 
picture, excitations that can be classified in accord with the centre group provide 
the disorder required to produce an area law of the Wegner- Wilson loop and, 
therefore, confinement. One striking feature is that — unlike monopole currents 
- centre vortices form gauge invariant two-dimensional objects, such that in four 
space-time dimensions, a linking number between a Wegner- Wilson loop and cen- 
tre vortices can unambiguously be defined, providing a geometric interpretation 
of the confinement mechanism J7|. 

I will restrict myself to discussing the superconductor picture which is based on 
the concept of electro-magnetic duality after an Abelian gauge projection and has 
originally been proposed by 't Hooft and Mandelstam ||. The QCD vacuum is 
thought to behave analogously to an electrodynamic superconductor but with the 
roles of electric and magnetic fields being interchanged: a condensate of magnetic 
monopoles expels electric fields from the vacuum. If one now puts electric charge 
and anti-charge into this medium, the electric flux that forms between them will 
be squeezed into a thin, eventually string-like, Abrikosov-Nielsen-Oleson vortex 
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Figure 1: Electric field distribution between two static SU(2) sources in the MA 
projection. Everything is in lattice units a « 0.081 fm. The sources are located 
at the coordinates (—7.5,0) and (7.5,0). 



which results in linear confinement. 

In all quantum field theories in which confinement has been proven, namely 
in compact U(l) gauge theory, the Georgi-Glashow model and SUSY Yang-Mills 
theories, this scenario is indeed realised. However, before one can apply this 
simple picture to QCD or SU(N) chromodynamics one has to identify the relevant 
dynamical variables: it is not straight forward to generalise the electro-magnetic 
duality of a U(l) gauge theory to SU(N) where gluons carry colour charges. How 
can one define electric fields and dual fields in a gauge invariant way? 

In the Georgi-Glashow model, the SO (3) gauge symmetry is broken down to 
a residual U(l) symmetry as the vacuum expectation value of the Higgs field be- 
comes finite. It is currently unknown whether QCD provides a similar mechanism 
and various reductions of the SU(N) symmetry have been conjectured. In this 
spirit, it has been proposed Q to identify the monopoles in a U(l) N ~ l Cartan sub- 
group of SU(N) gauge theory after gauge fixing with respect to the off-diagonal 
S'f/(iV)/?7(l) iV_1 degrees of freedom. After such an Abelian gauge fixing QCD 
can be regarded as a theory of interacting photons, monopoles and matter fields 
(i.e. off-diagonal gluons and quarks). One might assume that the off-diagonal 
gluons do not affect long range interactions. This conjecture is known as Abelian 
dominance ||. Abelian as well as monopole dominance are qualitatively realised 
in Lattice studies of SU{2) gauge theory || in maximally Abelian (MA) gauge 
projection, which appears to be the most suitable gauge fixing condition. 

In Figure [l], I display the electric field distribution between SU(2) quarks, 
separated by a distance r = 15a « 1.2 fm, that has been obtained within the MA 
gauge projection. Everything is measured in lattice units a « 0.081 fm where 
the physical scale derived from the value a/k = 440 MeV for the string tension is 
intended to serve as a guide to what one might expect in "real" QCD. Indeed an 
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elongated Abrikosov-Nielsen-Oleson vortex forms between the charges. In Fig. 0, 
I display a cross section through the centre plane of this vortex. While the elec- 
tric field strength decreases with the distance from the core, the modulus of the 
dual Ginzburg-Landau (GL) wave function, /, i.e. the density of superconducting 
magnetic monopoles decreases towards the centre of the vortex where supercon- 
ductivity breaks down. In this study the values A = 0.15(2) fm and £ = 0.25(3) fm 
have been obtained for penetration depth and GL coherence length, respectively. 
The ratio A/£ = 0.59(13) < 1/V2 corresponds to a (weak) type I superconductor, 
i.e. QCD flux tubes appear to attract each other. For details I refer the reader 
to Ref. Q 



3 String breaking 

In the pure gauge theory results presented above, the energy stored in the vortex 
between charges increases in proportion to their distance ad infinitum. In full 
QCD with sea quarks, however, the string will break into two parts as soon as 
this energy exceeds the energy required to create a quark-antiquark pair from 
the vacuum: inter-quark forces at large separation will be completely screened 
by sea quarks and excited T states can decay into a BB meson pair. In Fig. [3], 
I display a recent comparison between the quenched and nf = 2 static potential 
by the T\L collaboration at a sea quark mass m u d ~ m s /3 |[L1|| . Estimates of 
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masses of pairs of static-light mesons into which the static heavy-heavy systems 
can decay are also included into the figure. The potentials have been matched to 
each other at a distance r = ro ~ 0.5 fm. In presence of sea quarks anti-screening 
is weakened and, therefore, starting from the same infra red value, the effective 
QCD coupling runs slower towards the a s = ultra violet limit. This effect 
explains why at small r the unquenched data points are somewhat below their 
quenched counterparts: the effective Coulomb force remains stronger. Around 
r = 1.2 fm, the un-quenched potential is expected to become unstable. However, 
the data are not yet precise enough to resolve this effect. 

Motivated by such QCD simulations, the dynamics of string breaking has 
recently been analysed in some toy models First results on interactions 

between two B mesons in quenched QCD have been reported by Pennanen and 
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4 Glueballs and quark-gluon hybrid states 

In Fig. not only the ground state potential but also a so-called hybrid potential 
is displayed in which the gluonic component contributes to the angular momen- 
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Figure 4: The quenched glueballs spectrum (from Ref. p7| ). 



turn. Recently, the spectrum of such potentials has been accurately determined 
by Juge, Kuti and Morningstar |14|. The presence of gluons in bound states 
should also affect light meson and baryon spectra: one would expect additional 
excitations that cannot be classified in accord with the naive constituent quark 
model. On the Lattice and in experiment it should be most easy to discriminate 
states with exotic, i.e. quark model forbidden, quantum numbers from "standard" 
hadrons. Spin-exotic baryons cannot be constructed but only mesons and glue- 
balls. First results on light hybrid mesons have been reported by two groups [[L5[ . 
The lightest spin-exotic particle has quantum numbers J PC = 1~ + and a mass 
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between 1.8 and 2 GeV. Recent investigations incorporating sea quarks [16 
firm these findings. However, at present all experimental candidates have masses 
smaller than 1.6 GeV. Therefore, in the interpretation of experiment mixing be- 
tween spin exotic mesons and four-quark molecules, such as a 7r/i, should be 
considered. It is certainly worthwhile to investigate this possibility on the Lat- 
tice too. 

Recent results by Morningstar and Peardon ]T7J have revolutionised our knowl- 
edge on the quenched glueball spectrum. Only in case of the scalar glueball they 
fail to reach the precision of the 1993 state-of-the-art Lattice predictions [|TJ|]: 
finer lattices are required for a safer continuum limit extrapolation. As can be 
seen from Fig. 0, the ordering of glueball states has become fairly well established. 
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The fact that the lightest spin-exotic state 2 + ~ lies well above 4 GeV explains 
why such states have escaped observation so far. Glueballs are quite heavy and 
spatially rather extended, due to the lack of quarks that tie the flux together. 
Therefore, these states lend themselves to the use of anisotropic lattices: the size 
of the temporal lattice spacing is dictated by the heavy mass that one wishes 
to resolve while a much coarser spatial spacing can be adapted to resolve the 
glueball wave function. Introducing this anisotropy was vital for the improve- 
ment achieved. Recent results in QCD with sea quarks on the scalar and tensor 
glueballs are consistent with quenched findings Jl9|, [ll] . Beyond the quenched ap- 
proximation, glueballs will mix with standard quark model states. Investigations 



of such mixing and decay rates of the mixed states are challenging questions |g0 
that are waiting to be approached by Lattice studies in the near future. 



5 Light hadrons 

In addition to quantities that theorists or experimentalists are interested in, well- 
known observables can be computed on the Lattice too. The motivation is two- 
fold: testing QCD and gauging the Lattice methodology. Experimental low en- 
ergy input like the hadron spectrum is required in the first place to fix the lattice 
spacing and quark masses. Subsequently, among other predictions, the funda- 
mental parameters a s and quark masses || can be converted to, for instance, the 
MS scheme that is convenient for perturbative continuum calculations. It is not 
a priori clear whether the low energy results are compatible with values required 
to explain high energy QCD phenomenology. 

Assuming that QCD is the right theory, the observed states can serve as a 
guideline to judge the viability of approximations, such as quenching in the ab- 
sence of high precision full QCD results. Last but not least, quark masses and 
other parameters can be varied and Lattice results can be confronted with predic- 
tions of, for instance, chiral perturbation theory. Indeed, evidence for quenched 
chiral logarithms has been reported |21| . 



In Fig. [5], I display results from a recent state-of-the-art calculation of the 



quenched light hadron spectrum by the CP-PACS collaboration [21]. For com- 



parison the results from the GF11 collaboration |22| that have set the standard 
back in 1994 are included (squares). The 7r and p masses that have been used 
as input values for the lattice spacing a and quark mass m u = are omitted 
from the plot. m s has been set by two methods: forcing the K mass to agree 
with experiment (full circles) and forcing the mass to agree (open circles). Nei- 
ther of the methods can bring the spectrum completely in line with experiment. 
However, no mass comes out to be wrong by more than 10 %, indicating that 
the main effect of sea quarks is to renormalise the over-all value of the coupling, 
rather than altering mass ratios, despite the fact that all particles displayed, with 
the exception of the nucleon, become unstable in full QCD. First un-quenched 
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The quenched light hadron spectrum (from Ref. [21]) 



results by the same collaboration indicate an improvement in the direction of 
the experimental values. Many groups are at present studying quantities which 
one might expect to be more sensitive towards quenching like the rj mass, quark 



masses and the irNa term [p3 . 
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